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Abstract 



We have investigated the effect of friction on the DCC domain formation. 
We solve the Newton equation of motion for the 0(4) fields, with quenched 
initial condition. The initial fields are randomly distributed in a Gaussian 
form. In one dimensional expansion, on the average, large DCC domains 
can not be formed. However, in some particular orbits, large instabilities 
may occur. This possibility also greatly diminishes with the introduction of 
friction. But, if the friction is large, the system may be overdamped and then, 
there is a possibility of large DCC domain formation in some events. 
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The possibility of forming disoriented chiral condensate (DCC) in relativistic heavy ion 
collisions has generated considerable research activities in recent years. The idea was first 
proposed by Rajagopal and Wilczek They argued that for a second order chiral phase 

transition, the chiral condensate can become temporarily disoriented in the nonequilibrium 
conditions encountered in heavy ion collisions. As the temperature drops below T c , the chiral 
symmetry begins to break by developing domains in which the chiral field is misaligned 
from its true vacuum value. The misaligned condensate has the same quark content and 
quantum numbers as do pions and essentially constitute a classical pion field. The system 
will finally relaxes to the true vacuum and in the process can emit coherent pions. Since the 
disoriented domains have well defined isospin orientation, the associated pions can exhibit 
novel centauro-like ||-§| fluctuations of neutral and charged pions |PHH| . 



Most dynamical studies of DCC have been based on the linear sigma model, in which 
the chiral degrees of freedom are described by the real 0(4) field $ = (a, II), having the 
equation of motion, 

[□ + A(0 2 - v 2 )}$ = Hn a (1) 

The parameters of the model can be fixed by specifying the pion decay constant, f n =Q2 
MeV and the meson masses, m 7r =135 MeV and m cr =600 MeV, leading to A = (ml — 
m2)/2/2=20.14 and v = \(m 2 a -3ml) /(ml -m 2 )] 1 / 2 f n =&§.7l MeV and H = (120.55MeV0 3 
|T3| . It is apparent from eq.|l] that the vacuum is aligned in the a direction $ mc = (f^, 0) 
and at low temperature the fluctuations represent nearly free a and 7r mesons. At very high 
temperature well above v , the field fluctuations are centered near zero and approximate 0(4) 
symmetry prevails. 

It is instructive to decompose the chiral field, 

$(r,t) =< 0(t) > +S(f>(r,t) (2) 

where < <p >is the average over a suitable region of space and can be identified with the 
(local) order parameter and 5(f) are the semiclassical fluctuations and can be identified with 
quasi-particle excitations. Using eq.|2| and taking the average of eq.|l|, the equation of motion 
for the mean fields in the Hartree approximation can be obtained as, 

° 2 ^ f > = X(v 2 - < <j) > 2 -3 < 5(f)\ > - < 5(f) 2 ± >) < (j) > +Hn a (3) 

where 5<p» is the component of the fluctuation parallel to < <fi > and 5<f)± is the orthogonal 
component. This equation imply that the motion of the mean field is determined by the 
effective potential, 



V(< >) = ^(< > 2 +3 < 5<P\ > + < 5(f)l > -v 2 ) (4) 



which clearly differs from the zero temperature one in presence of fluctuations. By varying 
the fluctuations, chiral symmetry can be restored or spontaneously broken. It is also evident 
that the evolution of the order parameter critically depends on the initial values of the 
fluctuations. When 5 2 = (3 < 5(f) 2 > — < 5(f) 2 ± >)/6 is large enough the chiral symmetry is 
approximately (as 0) restored. If the explicit chiral symmetry breaking term is neglected, 
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the phase transition takes place at the critical fluctuations 5 2 = v 2 /6. For 5 2 < 8 2 , the 
effective potential takes its minimum value at < <fi >= (cr e ,0), where a e depends on 5 2 . 
When the mean fields are displaced from this equilibrium point to the central lump of the 
Mexican hat (< <fi >~ 0), the effective mass square 

m 2 = X(v 2 - < <f) > 2 -3 < 5<P\ > - < 5<f)l >) (5) 

will become negative and DCC can form. Since the domain size is directly related to the time 
scale, during which the effective mass remains negative, it strongly depends on the initial 
condition, 5 2 and < > of the system. In the one loop effective theory, the fluctuations can 
be replaced by (T 2 /2), their counterpart in a finite temperature field theory |[l^Jl5 . 



Recently, Biro and Greiner [0, using the Langevin equation for the linear sigma model, 
investigated the interplay of friction and white noise on the evolution of the order parameter. 
In the </> 4 theory, it can be shown explicitly that if the hard modes are integrated out on the 
two loop level, Langevin type of equation emerges for the propagation of long wavelength 
fields They found that in different realistic initial volumes ranging from 1 to 1000 

fm 3 , the average evolution do not shows any sensible instability. However individual events 
sometime shows significant growth of fluctuations. 

In the present paper, we will investigate the role of friction on the possible DCC domain 
formation. To this end, we solve eq.|], augmented by a Bjorken type of scaling expansion 



and Rayleigh dissipation term ||T9 

d 2 < 6 > ,D N <9 < 6 > , r9 , 9 „ , . 

dJ + (- + rj) = A[i; 2 - < > 2 -T 2 /2] < <\> > +Hn a (6) 

In the weak coupling limit, the friction r\ is related to the on-shell plasmon damping 
rate, r\ = 2^ v \. In the standard 4 theory, the plasmon damping rate can be calculated for 



the a and the II fields [0,|TH]. Assuming that the meson masses are the same, the friction 



coefficient can be obtained as [16|, 



V = llvi = 7#iA 2 — f Sp (l - e~ m ' T ) (7) 



where f Sp = - JT dtf\ is the Spence function. At T = T c = y/2f 2 - 2m2/A=123 MeV, 
and if m/T ~ 1, the friction 77 = 2.2fm~ 1 . This friction may be too high. Therefore we 
also investigate scenarios with 1/2 and 1/4 of its value. We also neglect the temperature 
dependence of the friction coefficient. 

We will solve the eq.|^ with quenched initial condition, for one dimensional (D=l) ex- 
pansion. As argued before, DCC formation depend critically on the initial condition. To 
reflect the uncertainty in the initial condition, the initial fields are randomly distributed to 
a Gaussian form with the following parameters, 

< <p >= (8) 

< (p 2 > - < <p > 2 = v 2 /6 (9) 

< >= (10) 

< I 2 > - < I > 2 = Av 2 /Q (11) 
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For each trajectories, we compute the characteristic quantity m t = sgn{m 2 )^\m\. The 
phenomenon of long wavelength DCC amplification will occur whenever this characteristic 
quantity is negative. We also calculate the amplification coefficient defined as, 

G = J \m\Q(-m 2 )dt (12) 

Amplification of unstable orbits are ~ exp(G). 

In fig.l, we have shown the ensemble average of 1000 trajectories of the characteristic 
quantity m t (solid line). We have considered four different values for the friction coeffi- 
cient, 77=0,0.5,1.1 and 2.2 fm~ 1 . Also shown is the m t for the most unstable orbit i.e. the 
orbit for which the amplification coefficient G is maximum (the dashed line). For 77=0, 
< m t > execute oscillation about its relaxed value of 140 MeV. It remain positive through- 
out the evolution. This indicate that in one dimensional expansion, on the average DCC like 
phenomena will not occur. The most unstable orbit shows nice oscillations with slowly di- 
minishing amplitude, indicating possibility of large DCC domain formation in this particular 
event. Thus while on the average, one dimensional system do not exhibit DCC phenomena, 
in some particular event, large DCC formation can occur. These results are in agreement 



with the other calculations p0l|T6f . With the introduction of friction the results changes. For 
small friction r/=0.5 fm^ 1 , < m t > do not oscillate much. It also remain positive through- 
out the evolution. Like the previous frictionless case, DCC phenomena is not expected on 
the average level. The most unstable orbit exhibit rapidly damping oscillation. It remains 
negative but for a short period of time. As size of DCC domain is directly related to the 
time during which mt remain negative, the result indicate that even for the most unstable 
orbit, large DCC domains can not be expected in this situation. Introduction of friction 
greatly diminishes the possibility of large DCC domain formation. With further increase of 
the friction result changes again. With 77=1.1 fm^ 1 , < m t > do not show oscillation. But 
we find that it enters into the unstable region for a short duration of time. With moderately 
large friction, DCC phenomena can occur on average level also. The most unstable orbit 
shows a interesting behavior. Once m t reaches a negative value, it remains negative for a 
long duration, then relaxing to its stable value of 140 MeV. Similar behavior is seen for 
77 = 2.2/m -1 . m t for the most unstable orbit stays negative for quite a long time (~15 
fm), indicating possibility of a large DCC formation in this particular event. The results 
indicate that if the friction is sufficiently large, one can hope to form large DCC domains 
in some particular event. This result can be understood easily. With large 77, the system is 
overdamped and once the system enters into unstable regime, it can not come out from it 
and remains unstable for long duration. 

In order to check that the oscillations in mt are stable or not, we have performed the 
time series analysis according to Hurst [0,0 an d obtained Hurst coefficient (H). Hurst 
coefficient varies between 0-1. For a perfectly random data set it is 0.5. If H > 0.5, the 
data set is persistent and they are antipersistent if H < .5. For the persistent data set, if 
the trend or behaviour in the data set is increasing or decreasing over a certain unit interval 
of time, it will continue to increase or decrease over such an interval, but no such comment 
can be made about antipersistent orbits. The most unstable orbits shown in fig.l are found 
to be persistent (H > .5). In fig. 2, we have shown the distribution of the Hurst coefficient 
for the 1000 trajectories calculated. For the frictionless system, quite a large number of 
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trajectories are found to be antipersistent. The oscillation for these orbits may eventually 
die out. With increasing friction, this number decreases and with 77=2.2 we find all 

the trajectories are persistent. 

The amplification coefficient (G) is an important parameter in the theoretical analysis of 
DCC. Larger the G, larger is the DCC domain size. In fig. 3, we have shown the distribution 
of G for all the 1000 trajectories as obtained in the present calculation. The distribution of 
the persistent trajectories (H > .5), are also shown (the shaded region). It is evident that 
maximum unstability is obtained in the frictionless case, with G as large as 2.4. Frictionless 
case also shows marked difference when only the persistent trajectories are considered. For 
the persistent trajectories the distribution decreases continuously with G, in contrast, when 
all the trajectories are considered, the number of trajectories remains nearly constant be- 
tween G=.5-1.5. With friction, as most of the trajectories are persistent, there is not much 
difference between the total distribution and the distribution of persistent trajectories. As 
is the case with frictionless case, we find that number of trajectories decreases continuously 
with increasing G. 

The pions from a DCC domain shows anomalous behavior, as they have definite isospin 
orientation. Defining the neutral to total pion ratio / 

7T° 

/= + ^ IT ( 13 ) 

the probability to obtain a particular fraction / in case of DCC formation can be obtained 

as 



Pif) = (M) 
incontrast to the binomial distribution peaked at 1/3 for normal hadronic reaction. From 



eq.|J, an equivalent probability that f < fi can be obtained, 

P(f <h) = ff (is) 

Noting that the pion yield will be proportional to / U 2 tcIt, we have calculated this 
probability. The results are shown in fig. 4 (solid line). The dashed lines are the probability 
distribution given by eq.fTH While for frictionless and for small friction, the agreement 



between the two are reasonable, for moderate to high friction, they in very good agreement. 
It seems that the pions are all from DCC like events. This raises a question. Preceding 
analysis shows that it is not possible to obtain large DCC domain except in a few trajectories. 
Yet we obtain neutral to total pion ratio following the probability distribution of a DCC 
domain. This indicate that 1 / v7 may not be good signal to identify the pions from DCC 
domains. 

To conclude, we have investigated the effect of friction on the possible DCC domain 
formation. It was shown that while in one dimensional case, on average DCC formation 
is not possible, there is finite probability for DCC domain formation in a event by event 
basis. We also find that friction strongly inhibit large DCC formation even in a event by 
event basis. However, if the friction coefficient is sufficiently large, the system may become 
overdamped and then large DCC domain can be formed. 
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FIGURES 



FIG. 1. Evolution of the ensemble averaged mt is shown for different friction coefficients (the 
solid line). The most unstable orbit, for which the amplification coefficient is the largest is also 
shown (the dashed line). 

FIG. 2. Histogram showing the distribution of Hurst coefficient for the 1000 trajectories. 

FIG. 3. Distribution of the amplification factor G. The shaded distribution is for the persisten 
trajectories with Hurst coefficient H > 0.5. 

FIG. 4. The probability that the neutral to total pion ratio / is less than f\. The solid line is 



the calculated while the dashed line is from eq,|l5 
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FIGURES 
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FIG. 1. Evolution of the ensemble averaged m t is shown for different friction coefficients (the 
solid line). The most unstable orbit, for which the amplification coefficient is the largest is also 
shown (the dashed line). 
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FIG. 2. Histogram showing the distribution of Hurst coefficient for the 1000 trajectories. 
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FIG. 3. Distribution of the amplification factor G. The shaded distribution is for the persisten 
trajectories with Hurst coefficient H > 0.5. 
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FIG. 4. The probability that the neutral to total pion ratio / is less than /]. The solid line is 
the calculated while the dashed line is for \pf probability 
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